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We study the properties of a spin-polarized Fermi gas in a harmonic trap, using the semiclassical (Thomas- 
Fermi) approximation. Universal forms for the spatial and momentum distributions are calculated, and the results 
compared with the corresponding properties of a dilute Bose gas. 
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I. INTRODUCTION 

Trapped degenerate atomic gases provide exciting op- 
portunities for the manipulation and quantitative study 
of quantum statistical effects, such as the strikingly di- 
rect observation of Bose-Einstein condensation. |]l]-g[ Al- 
though perhaps not as dramatic as the discontinuities 
associated with bosons, the behavior of trapped Fermi 
gases also merits attention, both as a degenerate quan- 
tum system in its own right and as a possible precur- 
sor to a paired Fermi condensate at lower temperatures. 
Q The ideal Fermi gas is an old and well-understood 
problem; there are many familiar systems where the non- 
interacting Fermi gas is a good zeroth-order approxima- 
tion. Unlike electrons in atoms and metals and nucle- 
ons in nuclei, however, atomic gases interact by predom- 
inantly short-range interactions whose effects are weak in 
the dilute limit. 

Trapped atomic gases provide an excellent laboratory 
for studying quantum systems in a controlled, confined 
setting. The quadratic potential of harmonic traps pro- 
vide a particularly simple realization of the confined 
Fermi system. At low temperatures, both the Bose ||] 
and Fermi Q gases are expanded relative to a classical 
gas at the same temperature; for fermions, however, this 
effect is due to the Pauli exclusion principle rather than 
atom-atom interactions. While in the Bose case a phase 
transition separates the degenerate and classical regimes, 
a trapped Fermi gas undergoes a gradual crossover be- 
tween the classical limit and the compact Fermi sea. 

We calculate the spatial and momentum distributions, 
the chemical potential, and the specific heat of a trapped 
spin-polarized Fermi g function of temperature. 

The properties of gases with two-or-more trapped spin 
states will be discussed elsewhere. We find that the 
properties of harmonically trapped gases with different 
particle numbers can all be described by the same uni- 
versal functions, after suitable scaling of variables. As we 
will see below, observation of the spatial distribution of 
the trapped cloud would provide an explicit visualization 
of a real-space "Fermi sea." 



II. DENSITY OF STATES 

We consider N spin-polarized fermions of mass M 
moving in an azimuthally symmetric harmonic potential, 
with a single-particle Hamiltonian 



Muji 



(1) 



where LOr and Wz = Aw^ are the trap frequencies in 
the radial and axial directions, respectively. Since the 
Pauli exclusion principle prohibits close-approach of spin- 
polarized fermions, we may to a first approximation ne- 
glect short-ranged interactions between them (see below) . 
The single-particle levels of eq. (0) are familiar: 
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where n^, Uy, and are non- negative integers, and the 
zero-point energy has been suppressed. For temperatures 
greater than the level spacing (/csT 3> huj), we may treat 
this discrete single-particle spectrum as a continuum with 
density of states 
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III. ENERGY AND LENGTH SCALES 

The chemical potential /x(T, N) is given implicitly by 
the equation 



N 



9{E) dE 



(4) 



At zero temperature the Fermi-Dirac occupation factor is 
unity for energies less than the Fermi energy Ep = fj,(T = 
0, N), and zero otherwise. A straightforward integration 
of eq. (^) then gives ^ 

Ep - nujr[Q\N]^^^, (5) 
which sets the characteristic energy of the atomic cloud. 
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The characteristic size g of the trapped degenerate 
Fermi gas is given by the excursion of a classical 
particle with total energy Ep in the trap potential: 



Rf = [2Ef/Muj:^] 
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where ar = {fi/MuJrY^''^ the radial width of the Gaus- 
sian ground state of the trap. For large A^, the width 
of the degenerate Fermi cloud is much greater than the 
quantum length Ur , and the Fermi energy is much greater 
than the level spacing of the trap, due to the Pauli ex- 
clusion induced "repulsion" between fermions. 

Similarly, we may define a characteristic wavenumber, 
Kp, which is determined by the momentum of a free 
particle of energy Ep'- 



Kp = [IMEp/h^]^/^ = (48iVA)i/^ cr-i 
= {A8NX/R%y/^. 



(7) 

(8) 



From eq. (||) we see that Kp is roughly the reciprocal of 
the typical interparticle spacing in the gas. 

As an example, consider spin-polarized ^Li. The radial 
frequency of this fermionic isotope of lithium in the TOP 
trap of ref. |l| would be Ur = 3800 sec^^; the character- 
istic ground state length ar is 1.6 /im. The TOP trap 
has an intrinsic axial/radial ratio A = VS. For N — 10^ 
atoms, the radius Rp is 25 /im; the typical interparticle 
spacing 1/ Kp is 0.1 /im. The Fermi temperature for this 
gas would be 3.5 /iK. 



IV. CHEMICAL POTENTIAL AND SPECIFIC 
HEAT VS. TEMPERATURE 

For general temperature, the chemical potential jj. must 
be determined numerically using eq. (||). We can find 
analytic expressions, however, in the limits of high and 
low temperature. For low temperature (ksT <C Ep) the 
chemical potential is given by the Sommerfeld expansion: 



(9) 



The third and higher-order terms in the Sommerfeld se- 
ries vanish since the density-of-states is a quadratic func- 
tion of energy. At high temperatures {i.e., in the classical 
limit ksT Ep), we find 



fi{T, N) = -fcBrin[6( 



fcsT 3 
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Numerical results for /i(T) / Ep are compared with 
these two limiting forms in figure ^ Evidently the 
low temperature approximation is quantitatively accu- 
rate below ksT/Ep ^ 0.5, and the classical expression 
holds above ksT/Ep - 0.6. 
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Low T Approx. 
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LU 




1 .00 



FIG. 1. The chemical potential DS. temperature. Both axes 
are scaled by the Fermi energy, which results in a universal 
curve that applies to all harmonically trapped Fermi gases. 

In eqs. (^) and (|o|), the particle number N enters 
IJl{T,N) only through the Fermi energy Ep. This result 
holds generally for all temperatures, as can be seen by 
casting eq. (Q) in dimensionless form by scaling E, fi, and 
1//3 by Ep. (The same conclusion holds for any density 
of states of the form g{E) = AE^ , for constant A, b.) 
Figure |l| is therefore a universal curve in the sense that 
it applies to harmonically trapped Fermi gases containing 
any number of particles. 

The specific heat ||l^ per particle of the trapped Fermi 
gas, 1/A^ 9C//9T|jv, is shown in figure |^. It is a mono- 
tonic function of temperature. For low temperatures the 
specific heat per particle is Tr'^kBiksT/ Ep); at high tem- 
peratures C — "iNks, as we approach the equipartition 
limit for particles in a harmonic potential. 
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FIG. 2. The heat capacity vs. temperature. The heat ca- 
pacity is scaled by ksN and the temperature by Ep. The 
classical result is shown by the dotted line. 
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V. SEMICLASSICAL (THOMAS-FERMI) 
APPROXIMATION 

Since the exact eigenstates of the harmonic potential 
are well-known, the properties of a harmonically trapped 
ideal gas can in principle be found directly by summing 
over these states. It is useful, however, to have approxi- 
mate forms for various observables that can be computed 
directly in the large N limit, where the exact sums be- 
come unwieldy. 

In the "semiclassical" or Thomas-Fermi approx- 
imation, [0 the state of each atom is labeled by a po- 
sition r and a wavevector k, which can be viewed as the 
centers of a wavepacket state. The energy of the parti- 
cle is simply the corresponding value of the Hamiltonian; 
the density of states in the six-dimensional phase space 
(r,k) is {2tt)^^, where sums over states are replaced by 
integrals over phase space. These semiclassical approxi- 
mations are valid in the limit of large N, as discussed in 
the Appendix. 

In the semiclassical limit, the number density in phase 
space is 



p{r,k■,T,^l) 
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(11) 



The chemical potential is given implicitly by the require- 
ment 



N 



J d^r d^kp{r,k;T,fi). 



(12) 



It follows from the correspondence principle that the 
Thomas-Fermi calculation of fi(T, N) using eq. (|l2|) re- 
produces the exact result obtained from eq. (^; this is 
easily confirmed for the harmonic oscillator, since the two 
integrals are related by a simple change-of- variables. 

After computing Iji(T,N), it is straightforward to cal- 
culate the spatial and momentum distribution functions 



n{k;T) = J d^r p{r,k;T,fj,). 



(13) 
(14) 



VI. 



SPATIAL DISTRIBUTION AT ZERO 
TEMPERATURE 



At zero temperature, we may define a "local" Fermi 
wavenumber kp{r) by 



2M 



V{v)^Ep, 



(15) 



where T^(r) is the trap potential. The density n(r) is 
then simply the volume of the local Fermi sea in fc-space, 
multiplied by the density of states (27r)~'^: 



n{Y-T = 0) 



kp{vf 
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(16) 



Note that n(r) vanishes for p > Rp, where p is the effec- 
tive distance 



p=[x^+y^ + X^z^]^/\ 
Combining eqs. (|l|) and (|l6|), we find §,|l2| 



, „ N NX 8 - .q/n 



(17) 



(18) 



for p < Rp. The cloud encompasses an ellipsoid with 
diameter 2Rp in the xy plane, and diameter 2Rp/X along 
the z axis. This aspect ratio is the same as that of a 
classical gas in the same potential, since the Boltzmann 
distribution is proportional to exp[—Mijj'^p'^/2kBT]. 



VII. MOMENTUM DISTRIBUTION AT ZERO 
TEMPERATURE 

One way to characterize the state of a trapped gas is to 
allow a rapid adiabatic expansion and then measure the 
velocity distribution by time-of-flight spectroscopy. 0] 
For the degenerate Bose gas, the observed anisotropy of 
this velocity distribution is dramatic evidence for quan- 
tum statistical effects. The semiclassical momentum dis- 
tribution for a degenerate Fermi gas at zero-temperature 
is simply 



n{k;T = 0) = 



1 



(27r)3 



dh Q{kp(r) - |k|), (19) 



where &{kp{r) — |k|) is the unit step function. The inte- 
gral (|l^ is the real-space volume within which the local 
Fermi wavevector exceeds Ikl: 



n(k;T = 0) 



N 8 



J^p 
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(20) 



where the maximum occupied wavenumber Kp was de- 
fined in eq. (0). Note that kp{r = 0) = Kp. 

Despite the spatial anisotropy of the trap, the momen- 
tum distribution of the degenerate Fermi gas is isotropic. 
This isotropy is a general feature of trapped Fermi gases, 
independent of trap potential, [ p^ since from eq. ( [l9| ) we 
see that h(k) depends only on the magnitude of k. 

The spatial and momentum distributions (^, 20) both 
have the same functional form, because 7i is a quadratic 
function of both position and momentum. In this sense, 
the distribution (|l^ ) can be viewed as a Fermi sea in 
real-space. The anisotropy of n(r) is due to the unequal 
spring constants, while the isotropy of n{k) is due to 
the isotropy of mass. That is, pi, Py, and pi enter the 



Hamiltonian with the same coefficient, while x , y" 



and 



need not. 
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VIII. NUMERICAL RESULTS 



In the semiclassical approximation, the spatial and mo- 
mentum distributions are easily determined numerically 
for any temperature as described in section As with 
the chemical potential, an appropriate scaling of these 
two distributions yields a universal form for all harmon- 
ically trapped Fermi gases when plotted vs. the scaled 
variables o/Rf and |k|/if^-, respectively. 

Figure shows the scaled density vs. scaled distance 
for keT/EF of 0, 0.25, 0.5, 0.75, and 1. At low tem- 
peratures, the density is close to their zero temperature 
form, with a thin evaporated "atmosphere" of thickness 
~ RpiksT / Ef) surrounding a degenerate liquid "core." 

In the classical limit, the density approaches a Gaus- 
sian in /9, with a width given by the equipartition theo- 
rem: {p^) = 'SRp{kBT/ Ep). As shown in figure ||, this ac- 
curately describes the density distribution for ksT/EF = 
1. The evolution of the density profile from its low- 
temperature Fermi form (|l8| ) to the classical limit can be 
tracked by calculating the mean-square excursion (/O^), 
which is shown in figure ^ in the dimensionless form 
(p^) / Rp ksT/EF. This is again a universal curve 
for all harmonically trapped Fermi gases. 

As noted above for zero temperature, the momentum 
distribution at non-zero temperature has the same form 
as the spatial distribution, since momentum and position 
both enter the single particle Hamiltonian quadratically. 
Thus the scaled momentum distribution _fi'|,n(k)A/V vs. 
|k|/iC_F is also given by figure ^. Similarly, figure ^ illus- 
trates the scaled mean-square momentum {'k^)/Kp vs. 
|k|/ifi? as well. 




P/Rp 

FIG. 3. The universal spatial and momentum density dis- 
tributions for ksT/EF = (bold), 0.25, 0.5, 0.75, and 1.0. 
The classical result for UbT/Ef = 1 is shown as a dashed 
line. 
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FIG. 4. The mean-square variation of the size of the cloud 
interpolates between a low-temperature degenerate regime 
and a high-temperature limit that is well-described by the 
equipartition theorem. 



IX. PERTURBATIONS 

What happens if the potential is not perfectly har- 
monic? We may treat 5V{y) as a perturbation. Here we 
focus our attention on the T = case. From cq. (|l5|), a 
change in trap potential shifts the local Fermi wavenum- 
ber by 



5kF{r) 



M 



kkplj) 



[SEF-SVir)], 



(21) 



where dEp is the change in Fermi energy. From eq. (p^), 
the corresponding change in density is 



Snir)^^^[6Ep-SVir)], 



(22) 



where the Fermi energy is adjusted to make / (Pr 6n{r) 
vanish: 



SEp = 



/ (frSV{r)kF{r) 
J d'^rkp{r) 



(23) 



Interactions between particles (which are small for 
spin-aligned gases) can be treated as a perturbation 
6V{r) = Un{r), where [/ is a parameter related to the 
p-wave scattering length. 



X. COMPARISON WITH THE BOSE GAS 

The interacting Bose gases of refs. and are in the 
Thomas- Fermi regime [|l4| . Since the gases remain dilute, 
two-body scattering may be treated by a delta-function 
pseudopotential of strength U — Airh'^a/M, where a is 
the s-wave scattering length. When the dimensionless 
parameter UN/TuvXa]^ is large (as is appropriate for the 
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experiments of refs. Q and 
interacting Bose gas is [Q 



the density profile of the 



with maximum radius 



Rb 



(24) 



(25) 



Note that the characteristic radius scales more slowly 
with particle number for the Fermi gas [N'^/^) than for 
the interacting Bose gas (A^^/^); similarly, the Fermi en- 
ergy scales as N^/^ while the zero temperature chemical 
potential of the Bose gas varies more rapidly, as N^^^. 

The axial/radial aspect ratio for both classical and de- 
generate trapped gases is A, since in every case the den- 
sities are functions of p only. The velocity (momentum) 
distributions, however, can be quite different: for classi- 
cal and Fermi gases the velocity distribution is isotropic, 
while n(k) for a zero temperature Bose gas is the square 
of the Fourier transform of ^/^^(r). This is notably 
anisotropic. Note that as N increases, both Rp and Rb 
increase, but the widths of the respective momentum dis- 
tributions go in opposite directions: Kp increases with 
iV, while the typical momentum of a particle in a trapped 
Bose condensate decreases with particle number, since 
Kb ~ ^/Rb by the uncertainty principle. 

It is amusing to compare the interatomic repulsion in 
a Bose gas with the effective repulsion experienced by 
fermions due to the Pauli exclusion principle. Equat- 
ing the characteristic Fermi and Bose radii (^ |2^) we 
see that, crudely speaking, the spatial distribution of 
a degenerate Fermi gas is mimicked by that of a Bose 
gas interacting via an effective "Pauli pseudopotential" 
UcS Ef{Rp/N), which is the characteristic energy 
multiplied by the volume per particle. Equivalently, the 
effective scattering length brought about by the Pauli 
principle is a^s ~ Kp^, i.e., the inter-particle spacing. 
(This is quite expected, since the inter-particle spacing 
is the only appropriate length in the ideal Fermi gas.) 
The use of such an effective interaction is limited by the 
fact that (a) the momentum distributions of the Fermi 
and Bose gases remain quite different and (b) the gas is 
not dilute with respect to the exclusion-induced "inter- 
actions" since Kpace is of order unity. 



XII. APPENDIX: VALIDITY OF THE 
SEMICLASSICAL APPROXIMATION 

The semiclassical approximation can be safely applied 
to an inhomogeneous Fermi gas of density n(r) if we can 
imagine partitioning the system into cells of linear di- 
mension £ such that the following two conditions are si- 
multaneously met: 

(1) The number of particles in a cell is much greater 
than unity, so that a local Fermi sea may be envisioned: 



n{r)f > 1. 



(26) 



(2) The variation of the trap potential across the cell 
{£WV) must be small compared with the local Fermi en- 
ergy fi^kpirY /2M , so that within a cell the potential 
energy is nearly constant. At low temperature, this con- 
dition becomes 



2Ar ^ " 



(27) 



where we have used eg . (|T^). 

Combining eqs. (26) and (27), we see that to be able 
to choose a (possibly r-dependent) cell size £ that simul- 
taneously satisfies these two conditions, the number of 
particles per quantum volume must satisfy 



n(r)(T'^ » -, 



(28) 



where a is as before the quantum length {h/MujY^^ and 
we have omitted factors of order unity. 

At low temperatures, the semiclassical density given 
by eq. (|l|) scales as N/R% ~ N^/^/a^ near the ori- 
gin, so the Thomas-Fermi approximation is always self- 
consistent at the center of the trap for large N. (This 
can be confirmed at r = by direct summation of the 
squares of the simple harmonic oscillator eigenfunctions 
up to energy Ep.) Near the periphery of the cloud, how- 
ever, the density becomes small, and the approximation 
is not valid. It is easy to show that semiclassical treat- 
ment fails within a shell at the periphery of the cloud, 
whose thickness SR ~ l/Kp ^ aN^^^^ vanishes in the 
limit of large N. Within this shell only the exponential 
tails of a few single particle states contribute to the den- 
sity; this is analogous to the corresponding region of the 
Bose gas. 15 
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